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The existence of a classical limit describing interacting particles in a second-quantized theory of 
identical particles with bosonic symmetry is proved. This limit exists in addition to a previously 
established classical limit with a classical field behavior, showing that the limit h — > of the theory 
is not unique. An analogous result is valid for a free massive scalar field: two distinct classical limits 
are proved to exist, describing a system of particles or a classical field. The introduction of local 
operators in order to represent kinematical properties of interest is shown to break the permutation 
symmetry under some localizability conditions, allowing the study of individual particle properties. 
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Quantum Field Theory was born in the attempt to conciliate Quantum Mechanics and Relativity and in the attempt 
to deal with the so-called "particle- wave duality" underlying quantum phenomena by making the corpuscular character 
of matter compatible with the classical notion of field. Its main purpose is to describe particle physics, but the way 
in which its fundamental principles admit or recognize the notion of particle is rather indirect, if not obscure. A 
fundamental attempt to understand when quantum field theories describe particles was the work of Haag and Swieca 
[l|, further developed by Buchholz and Wichmann in [2j and in subsequent works. In a nutshell, these works point 
to the fact that any relativistic quantum field theory describing particles must have some specific limitations on the 
number of degrees of freedom at finite volume and limited energy. 
Oh' Quantum Field Theory is believed to be, in some sense, a fundamental theory, but the notions of particles and of 
fields are derived from our sensorial experiences in a classical macroscopic world. It is, therefore, of great importance 
to have a precise understanding on how the classical scenario can be reached from that more fundamental quantum 
starting point. In particular, one should naturally expect that the particle-field duality manifests itself in any general 
attempt to reach the classical limit of quantum fields. The existence of these two different limits (particles or fields) 
lies deeply in the structure of quantum field theory and its physical interpretation. It was first remarked in the 
fundamental work of Hepp Q on the classical limit of quantum systems and the main purpose of the present work 
is to clarify certain aspects of this remark, specially in the relativistic regime. We believe that the analysis of these 
^{■^ . multiple classical limits could have some conceptual importance in the context of quantum field theories formulated 
in curved spacetimes, where no natural concept of particle states is available. 

There are several ways to approach the formulation of a classical limit of Quantum Mechanics. In this work we 
follow the ideas introduced by Hepp in [3| , which can be applied to a wide range of systems and can be understood in a 
simple and precise way. The central result of his work combines an old observation of Schrodinger Q , which led to the 
discovery of coherent states d, @ , with the intuitive explicit content of Ehrenfest theorem. Schrodinger observed that 
in a harmonic oscillator a Gaussian wavefunction moves without distortion along a classical orbit, what led him to 
try to understand particles in general as localized wavepackets in motion. But that could not be accomplished, since 
wavepackets in general spread, so that particles cannot remain localized indefinitely. On the other hand, Ehrenfest 
theorem describes the time-dependence of quantum averages with classical-looking equations. The main drawback is 
that for a function of the position operator A(q) one has (A(q)) ^ A((q)), unless A is linear. In the case where A is 
an external force, if one considers (q) as representing the position of the particle, a force A((q)) would be expect to 
act on it, whereas by Ehrenfest theorem the force is (A(q)). Thus (q) deviates from the classical trajectory. 

In order to briefly describe Hepp's analysis, consider a simple example. The classical nonrelativistic motion of a 
single particle of mass m moving in one-dimension under an external potential V is described in Classical Mechanics 

by the Hamiltonian H = ^ + V(x), leading to the canonical equations p = 3 = - 2 -< Let (£(t),7r(t)) be the 

solution for initial conditions x(0) = £, p(0) — n. The dynamics of the corresponding quantum system is defined by 
the Hamilton operator 

h 2 d 2 
2m oq A 

where q is the position operator. Let U (t) be the propagator associated to this Hamiltonian. The question we face is 
how to recover the classical trajectory defined by £(t) and n(t) from the quantum Hamiltonian. Following Hepp, there 
are three steps involved in the solution of this problem. In the first step we go to the Weyl representation, replacing 
position and momentum operators by their exponentiated versions exp(iaq), exp(ifep), a, b S R, so that only bounded 
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operators are involved. In the second step we consider averages of general Weyl operators W(a, b) := Q l i a i+ b P) on 
time-evolved coherent states \a), with a = (£ + in)/(\/2h). Finally, the limit % — * is taken. As established in Q, it 
turns out that under natural regularity requirements, and for t restricted to a finite interval \t\ < T, one has 

lim (o| U(t)*W(a, b)U(t) \a) = e^W+M*)] , (2) 

from which the classical trajectory (£(£), 7r(i)) can be recovered (for the above mentioned time- interval). The result 
is valid for a large class of potentials, as stated in 0]. 

The steps and results leading to can be easily extended in order to indicate how the classical limit of a quantum 
nonrelativistic system with finitely many degrees of freedom describing distinguishable particles has to be performed 
in order that this limit describes a classical mechanical system of finitely many particles (at least for a short time) . 
In this case one can use coherent states as those of Eq. ^ in connection with individual position and momentum 
operators pi, qi for each particle to find the classical limit of interacting point particles. The procedure presented 
above have also been extended in [3j to some quantum systems with infinitely many degrees of freedom, leading to 
classical limits describing classical fields. 

An important quantum system considered in [3J] for which the classical limit is a classical field is a a system of second- 
quantized interacting bosonic particles moving in one dimension, described in a bosonic Fock space T S ^H) = ©^Lo^"^ > 

(n) 

where H s is the usual Hilbcrt space of symmetric square- intcgrablc functions over R™, with the Hamiltonian 

H = ~J dxa*(x)V 2 a(x) + ±[ dxdya*(x)a*(y)V{x - y)a(x)a(y) , (3) 
where the creation and annihilation operators a* and a satisfy the usual commutation rules, 

[ a ( x ), a (y)} = i a *( x ), a *(y)} = °> i a ( x ), a *(y)} = $(x-y) , 

and V(x) = V(—x) — V(x)* is some real Kato potential. The classical limit of this system was studied in Q as a 
model for superfluid Helium, and that work was one of the motivations for Hepp's. In this case the classical field 
behavior is taken as an approximate description for the dynamics of the matter density in the fluid. That limit can be 
achieved with the same strategy adopted for the case of the one-particle dynamics of |T]). Coherent states of bosonic 
systems are usually defined in order to make field aspects become more evident, as can be seen from many examples 
in Quantum Optics (see e.g. Q for an introduction). Accordingly, if one considers the time-evolution of the average 
of exponentiated smeared local fields evaluated on an initially coherent state, one gets a classical field theory in the 
limit h — > 0. As shown in 0], the classical field a satisfies the non-linear partial integro-differential equation 

dot i f 

7ft (ft*.*) = ^V 2 a(f3,t,x)+i J dyV{x-y)\a((3,t,y)\ 2 a((3,t,x), (4) 

where (3(x) is the initial condition (which must be taken in D(V 2 )), and the constant /i is related to the mass m of 
the bosonic particles. The particular case when V{x — y) = gS(x — y) leads in ([?]) to the well-known Gross-Pitaevskii 
equation (or non-linear Schrodinger equation), widely employed in the study of Bose-Einstein condensates. The field 
behavior exposed in Q is not expected to hold along an arbitrary time interval. 

Since ([3]) is assumed to describe a quantum system of interacting particles, it is natural to expect that a second kind 
of classical limit exists which describes classical particles instead of classical fields. In fact, Hepp observes that these 
limits should depend on the way in which certain physical parameters are scaled when h — ► 0. There are, however, 
some implicit difficulties in applying Hepp's program to systems of indistinguishable particles, as in the case of the 
nonrelativistic quantum many-body system described in ([3]) or a relativistic quantum field model. The main problems 
are that : (i) the coherent states are not invariant under permutation symmetry, and (ii) observables describing 
individual kinematical properties may not be defined. The problem (i) is circumvented with the use of symmetrized 
products of single-particle coherent states, but the second problem requires some variation of the general technique. 
We will show how the introduction of local operators acting on a class of essentially localized states defined herein can 
be used to address this problem. A classical limit of N interacting particles is found when TV apparatuses situated 
at distinct regions are used to observe an essentially localized state with N localization centers coinciding with the 
apparatuses' positions. 

In the relativistic regime, one has to face the additional problem of defining the notion of quantum particle-like 
states in a proper way, since this is required for a particle classical limit. We adopt the notion of essentially localized 
states discussed by Haag and Swieca in [l[ as a suitable representation of the intuitive idea of a particle, and adjust it 
to our purposes. Position averages will be evaluated with the Newton- Wigner operator Q. An explicit construction 
based on the single-particle relativistic coherent states of [lfj will be shown to lead to the desired particle classical 
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limit of a free massive scalar field. As an intermediate step, we discuss the classical limit of the first-quantized Klein- 
Gordon field and show that the methods of [|[ can be adapted to solve this problem. Since the field classical limit of 
this theory is proved in 0, we conclude again that the classical limit is not unique. 

II. COHERENT STATES AND ESSENTIALLY LOCALIZED STATES 

Coherent states \a) for a one-dimensional quantum mechanical system are usually defined as eigenstates of the 
annihilation operator, a\a) — a \a), with a = (£ + z7r)/\/2fo € C It follows from this definition that \a) is a minimum 
uncertainty state centered at x = £ and p = tt with equal uncertainties Ap = Aq, and this property is one of the main 
motivations for the study of such states, since it allows one to think of them as the "most classical states" in some 
natural sense. An alternative and equivalent definition is given in terms of the action of a translation operator U(a) 
on the harmonic oscillator ground state |0), 

| a) := U (a) |0) , U (a) := exp(aa* — ad) , (5) 

where one has U(a)aU(a)* = a — a. The scalar products are given by 

|(a|/?)| 2 = exp(-|a-/3| 2 ), (6) 

thus coherent states are not orthogonal; nevertheless, the overlap decreases rapidly with the distance \a — (3\. 

Coherent states can be defined in much more general systems (see [ill The case of n-dimensional systems is 

straightforward, as well as the case of many-particles systems when no symmetrization is required. Let d be the number 
of spatial dimensions, and n the number of particles. In this case the coordinate space is TV-dimensional, with A = nd. 
All that is needed is to define a series of A annihilation operators a r = (q r + ip r )/y2h, one for each spatial dimension 
and particle of the system, and consider the simultaneous eigenstates of all annihilation operators, a r \a\, . . . , ajq) — 
a r \a\, . . . , ckjv), t = 1, . . . , N. The eigenstates |ai, . . . , ajj) are the coherent states for this system. Writing the 
labels a in terms of real and imaginary parts as a r — (£ r + wr r )/\/2ft, one finds that \a,\, ajq) describes a 
Gaussian wavefunction centered at x = . . . , £jv) in coordinate space and p = (jn, . . . , ttn) in momentum space 
with minimum uncertainty for each canonical pair, Aq r — Ap r = \J h/2. So everything goes as in the previous case. 

The case of many-body systems of identical particles, as the one described by the Hamiltonian ([3]) or in Quantum 
Field Theory, demands special care. Consider an A-particle state of a many-body system of identical particles in one 
space dimension and with bosonic symmetry (to which we restrict ourselves in this work). Because of the permutation 
symmetry, one cannot introduce individual position and momentum operators for each of the particles as in the case 
of distinguishable particles and a different procedure is required. The most natural choice, and the one adopted in 
this work, is to consider symmetrized products of single-particle coherent states \aii), 

\ai, a N ) s = T?y^ I ^(1)1 a *(N)) = ^ I'M 1 )) ® ' " • ® | a 7r(JV)) > (7) 

7T 7T 

as the analog of A-particle coherent states, where the sum is over all elements of the permutation group of A elements 
and where the normalization constant is 

M = N\y] (a^a^i)) ■■■ (a N \a^ N -)) . (8) 

These states will be of special interest in what follows, as examples of what will be called "essentially localized states" . 
It will be argued that such states reproduce within nonrelativistic quantum mechanics some basic intuitive properties 
of states so denoted in Algebraic Quantum Field Theory [l|. 

In order to discuss localization properties consider the example of one-dimensional single-particle states first, with 
the Hilbert space Ti = L 2 (R, dx). Denote by (A)$ := (^'l-d-IV') the expectation value of an operator A acting on TL in 
some vector state for tp = \ip) E TL with ||^|| 2 = (ip\ip) = 1. 

Let A be a bounded operator acting on Tt with norm ||A||. For any finite open region O C R, define the local 
version of A corresponding to the region O by 

Ao ~ \{xoA + A Xo ) , (9) 

where \o is the characteristic function of O, i.e., xo( x ) = 1 if a; G O and Xo(%) = if x ^ O. The operator Ao is 
also a bounded operator acting on TL and one has A — Aq — \ [(1 — Xo)A + A(l — xo)] ■ Hence, for expectation values 
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of A — Aq on some normalized vector state ip one has 

(A)^ - {Aoh = 5<(1 ~ Xo)1>\A\l>) + \{i>\A\{l - X o)V) , 
from which we get, by the Cauchy-Schwarz inequality, 

\{A)tp-{AoU\ < p||||(i- Xo )V>|| • 

Let us consider the concrete case where O = — R, £ + R), the interval of radius R > around the center £ S R. 
We say that ip e 7Y is essentially localized around £ if ||V'l|o c : = II (1 — xo)^!! j the fraction of the norm of the state 
lying outside O, satisfies 

U\\o* := ll(l-Xo)^|| < ICe-*''* 1 , (10) 
for all i? > 0, /C being some positive constant. For vector states satisfying (fTQ|) . one has 

((^-(Ao^l < IC\\A\\ c - R2 / 2h . (11) 

Thus the operator Ao can in fact be understood as a local version of the operator A in the region O, since the 
expectation (Ao)^ is a good approximation for (A)$, up to an error which decreases rapidly with R. 

In particular, this approximation works fine for coherent states. It was already mentioned that the position wave- 
function ip a (x) — (x\a), corresponding to a coherent state \a) is a normalized Gaussian curve with width y/h centered 
around some £ G R. Thus ip a {x) is essentially contained in a radius of a few \fh around £ and one has 

IWIo" : = ll(l-Xo)^a|| 2 = 1- / |Va(x)| 2 dai < Ke~ R (12) 

for some K, > 0, and decreases faster than exponentially with increasing R. 

There is another sense in which (TTU)) leads to a notion of essentially localized states. Consider two distinct states 
and \(f>) which satisfy (|10p with respect to disjoint regions O and O' , and let the distance between such regions 
be d(0,0') = d > 0. Hence, | (ip\<f>) | < 2e- d ^ 8R . Therefore, if and \cj>) are states essentially localized in disjoint 
regions, the overlap between the corresponding wavefunctions is small for large distances, and decreases as a Gaussian 
with the distance between the localization regions. Accordingly, the inequality (|10[) is used herefrom as the defining 
property of a "state essentially localized in O" . 

The concepts of local operator and essentially localized states can be easily translated to TV-particle states. The 
symmetrized product of single-particle states essentially localized in pairwise disjoint regions Oi (with OiDOj = 
if i j) is the natural extension adopted here. These states will be denoted "multiply localized states" or "essentially 
localized states" . Simple examples are given by the symmetrized product of coherent states as displayed in Eq. ([7]) , as 
long as the localization centers £j of the coherent states |oi) are chosen sufficiently far apart. Given a local operator 
Ao acting on the single-particle subspace we define its ./V-particle counterpart acting on Ti.i N ^ as 

N 

=5Z 1<8> •••1®A)«>1«> (13) 

i=l 

where the operator Ao occupies the z-th slot of the tensor product in the i-th term of the sum. 
The average value of the local operator A^ on a multiply localized state \ip) is given by 

( A o } ), = EE jf2 vWVV(i)) • • • mAo K W ) • • • (V'jvI^w) ■ (14) 

i 7T 

The factor N\/M 2 goes to 1 when K— * 0. In the case where \ip) is a multiply localized state with localization centers 
Oi, one has 

MVvw) < 2e- d2 / 8r \ if i^ir(i), (15) 

where d is the smallest distance between any two Oi, Oj, i ^ j. Thus, only a few terms in (|14| contribute appreciably 
to the sum, i.e., those for which i = Vi. Each such terms contains a factor (ipi\ Ao \tpi}- Now, suppose some of 
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the localization centers coincides with the region where the local operator in question is defined: Ok = O for some k. 
Then one has the bounds 

\{AU k -{AoU k \ < ||A||e-* 2 /\ (16) 
K^okl < \\A\\e- d2 / h , ifi^fc. (17) 

We assume in addition that 

\(A\A \^)\ < \\A\\e- d2 / 2h , ift^j. (18) 

This inequality corresponds to the requirement of a local property of the operators. It is automatically satisfied when 
Ao is the local version of any operator A of the form A — A(q). It is also satisfied when A = A(q) exp(ibp/H), with 
b < d, that is, when the operator A is a function of the coordinates, up to some small translation. In this case, we 
get from (fl4|) the approximation 

« {A) 4h , (19) 

with an error which decreases as exp(— d 2 /AK) exp(—R 2 /h) for large arguments d, R. Therefore, the average value of 
a local many-body operator associated with a region O, evaluated on a multiply localized state which has O as one 
of its localization centers, reduces to the average value of the corresponding single particle operator evaluated on a 
state essentially localized above O. 

The intuitive picture underlying the definitions and approximations given in this section is the following. In a many- 
body system of identical particles, one cannot in general disentangle kinematical properties of individual components, 
since permutation symmetry mixes them. But for some special states, an approximation is feasible where the particles 
behave as independent subsystems. Such "essentially localized states" look like isolated lumps of matter distributed 
over distinct regions of space. If a measuring apparatus works in a region where no lump of matter is present it will 
detect nothing. If there is one lump of matter present, a single particle will be detected, but the observer will be 
unaware of the existence of other identical particles comprising a larger Hilbert space together with the detected one 
- his "local operators" break the permutation symmetry since they are not sensible to particles far apart. If the 
measuring apparatus operates in a region where several lumps of matter are present, the individual properties of the 
particles become intertwined. The configurations of interest for the existence of a classical limit describing particles 
are those in which N observers situated at distinct regions can measure properties of one particle each. 



III. CLASSICAL LIMIT OF INTERACTING BOSONIC PARTICLES 

A vector state in Hf N can be represented by a normalized wavefunction ip(xi, . . . , xn) symmetric under the 
exchange of any two coordinates Xj, Xj. Let |a^) be single-particle coherent states and ipi(x) be the corresponding 
wavefunctions. The symmetrized product of such states is represented by 

tp(xi,...,x N ) = -^^ip^ixi) . . .tp^tf)(x N ) , (20) 

TT 

where J\f is a normalization constant. The action of the Hamiltonian Q on such states is given by 
h 2 1 

Hip(xi,...,x N ) = - — —^2^2\' 2 j ij 7T(1) (x 1 )...i(j 7liN) (x N ) 

j TT 

+ 2N r \'52^2 V ( Xj ~ Xk W< l )( Xl ) • ■ -VvwOsiv) • (21) 

Since the Hamiltonian is time-independent, the propagator is simply U(t) = exp{iHt/H). The potential V{x) is 
required to satisfy the regularity condition J da; ^(x)! 2 exp(— px 2 ) < oo for some p < oo [3(. This quantum dynamics 
is to be compared with the classical Hamiltonian equations 

H -Ey + ^E^-W, (22) 

(j=*i, rfi=-E^-&). (23) 
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Following [3J, we assumed that V 2 V is Lipschitz, so that solutions of the canonical equations exist and are unique 
in a finite time-interval \t\ < T. Denote by £(a,i), n(a,t) the solution for initial conditions oij = (£j + iiTj)/y/2H, 
3 = 1, N, and let V be Holder continuous C 2+c in an open neighborhood of £(a,i), for all \t\ < T. 

In order to compare the classical and quantum equations, consider for each 3 the localized operators W(a, &)c>j(t) 
associated to the Weyl operators W(a, b) := expi(ag + bp), acting on the single particle space, given as in (S]) by 



W(o, b) 0j {t) ■= \{xo } {t) expi(aq + bp) + expi(aq + bp)xo s (t) 

where Oj (t) := (£j (a,t) — R, £j (a, t) + R), where xo, ft) is the characteristic function of the set Oj (t) , and q, p are the 

and momentum operators. For each 3, let W(a, 6)^ 
as in (fT3|) . The desired classical limit is encoded in the expression 



usual position and momentum operators. For each 3, let W(a, b)^^ be the corresponding iV-body operator defined 



lim 

ft->0 



(i>\ U(t)* W{a, b) l ff (t) U{t) |V> - e^W+^W] 







(24) 



valid for each J = 1, . . . , N, where \ip) is the state defined in (|2"0)) . This limit is proved as follows. 

Consider the left side of the equation The state \tp) is a superposition of N\ unsymmetrized coherent-states 

in 7i® N . Let \<j>) be one of these states. The Hamiltonian (|2Tj) can be understood as acting on TC® N , and in this case 
it is known from Hepp's work Q that, under the assumed hypotheses, 



ft->o 



lim U(t) \(j>) - U(a(t))W{t) 



= 0. 



where W(t) = Tcxp —i/H J dt'H'(t') (here, T denotes the usual time-ordering prescription) is the propagator 
associated with the second-order Hamiltonian 



3 i¥=k 

obtained from the linearization of the operator in (|2ip around the classical orbit £(a, t), ir(a, t). The state |0) is 
the coherent state centered at zero position and momentum. To the state W(t) |0) there corresponds a wavefunction 
which is Gaussian and centered at zero in each coordinate Xi (an explicit expression is derived in Appendix [SJ . Thus 
one can write 



4>'{x, t) = (x\ U(a(t))W(t) |0) = Ylfiixj, t) 

3 



exp 



^u 3 {t){x 3 ~ &(a, t)) 2 + t ' x j 



(25) 
(26) 



The positive coefficients u>i(t) are continuous functions on [0, T] determined by the Hamiltonian H' . The state U(t) 
is the symmetric part of U(t) \<p), which for small h is well approximated by 



ip(x, t) 



1 

AH 



(27) 



with an error of order 0(h 5 ' 2 ). Define fi := min {wj(f)| 3 = 1, • ■ ■ , N and t £ [0, T]}. Then (|26p ensures that the 
states (|27p are essentially localized around the classical orbit £(a,i), i.e., ||<^(ie, £)llo c (t) < e ~ R ^ 2h , with h! = h/fl. 

Now observe that the Weyl operators can be rewritten as W(a, 6) := exp(iabh/2) exp(iaq) exp(i(bh)p j H) , so there 
is some fi such that bh is less than the smallest distance between the particles. Thus the inequality (I18|) is valid, and 
then the approximation (fT9|) yields, for each k = 1, . . . , N , 

(i>\ U(ty W(a, 6)W (t) U(t) IV-) « (0' fc (t)| W(o, b) \<f>' k (t)) ■ 

This approximation is valid as long as the classical orbits £j -(a, t),j = l, N, do not cross, i.e., \£j(a, t)—£i(a, t)\ > 
2R + k, for all j ^ I and t € [0, T], k > 0, in order that a finite minimal distance re exists between the localization 
regions centers )j(ct, t). The error involved in the approximation vanishes for h — ► 0, and the average values on the 
right side can be easily evaluated, leading to the limit (f2"4"|) , and completing the proof. 
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Some comments are in order. First, it is clear that there is a lot of freedom in the choice of the local operators, 
since there are many ways to define localization regions Oj(t). In the proof displayed above, these regions follow the 
classical trajectory with some fixed radius R. This is not necessary — it is sufficient that the classical trajectory of 
the j-th particle is in the interior of Oj(t). If there are disjoint regions Oj such that £(t) £ Oj for all t, then the time- 
dependence can be removed. In this case one would have a physical situation where the measuring apparatuses are 
placed at fixed regions where the particles are confined, a situation which is likely to happen for small time-intervals. 
Another remark concerns the possibility of collisions between particles. In this case the classical limit does not exist 
in the sense of the proved theorem, so the classical limit of scattering processes are not considered here. 

Regarding the existence of two distinct classical limits in the same quantum system, it is seen that the existence 
of two kinds of coherent states for bosonic systems is responsible for the fact. The usual field coherent states lead to 
a classical field equation obeying (0J for a time- interval of order t ~ ft 2 , while the essentially localized states defined 
herein lead to a particle dynamics for a time- interval of order h s / 2 . The introduction of local operators is necessary 
in the case of the particle limit in order to break the permutation symmetry and allow for the study of individual 
trajectories. 

A natural question connected with Hepp's analysis is to what extent it helps understanding the existence of a 
'classical world' as a consequence of more fundamental quantum laws. In short, what it states is that under a certain 
condition - the existence of a coherent state - and for a small time-interval, operator averages obey the expected 
classical laws of motion. Now, this result raises two natural questions. The first is to understand why a classical 
behavior is usually met with in macroscopic scales, that is, why states other than coherent states usually do not show 
up in macroscopic scales. Another question concerns the possibility of removing the restriction to small time-intervals, 
in order that this classical limit exists in time-scales compatible with all classical dynamics. 

A possible improvement of the theorem in order to address these problems could be the inclusion of an external 
system interacting with the system of interest. That would bring some contact with the widely studied decoherence 
approach to the emergence of classical behavior in quantum systems flij . It is known that several new effects can show 
up in this case, such as Zeno effect [l5| and non- unitary time-evolution [2], for instance, and it is not unreasonable 
that they may play an important role in the existence of a classical limit for large time-scales. Some experimental 
evidence points in this direction. A clean example is found in a series of papers on interference of matter waves 
(fullerenes) in a Talbot-Lau interferometer (UGH- ^ was observed that both interaction with gas particles as well 
as emission of radiation, i.e., interaction with a quantized electromagnetic field, helps preserving classical behavior. 
In both cases a process of localization of the particles is supposed to happen with a certain frequency in consequence 
of the interaction with the external system, thus naturally enforcing the regular occurrence of the conditions required 
for the validity of Hepp's analysis. 

IV. CLASSICAL PARTICLES AND THE KLEIN-GORDON FIELD 

In this section we extend our results to the relativistic regime. We will discuss the existence of two distinct kinds 
of classical limits in a system of bosonic relativistic particles, one of them describing a classical field, the other one 
describing classical systems with finitely many particles. For simplicity, we consider a scalar field theory with mass 
m > in 1 + 1-dimensional Minkowski spacetime. The existence of a classical field limit for this system in the presence 
of a polynomial interaction was proved in [H, roughly in the same way as for nonrelativistic mechanics: one studies 
the time-evolution of average values evaluated in an initially coherent state, and verifies that these averages obey 
the expected classical equations of motion. Here we show that a classical particle limit can also be reached, using 
methods analog to those applied in the nonrelativistic case, i.e., by studying the time-evolution of the average of local 
operators evaluated on essentially localized states. This will be done explicitly for the case of the relativistic position 
operator, in order to show how classical trajectories arise from the quantum dynamics. 

It is well-known that the problem of localizability of relativistic particles is rather more intricate than in the 
nonrelativistic case. The basic difference is that here the construction of a wavepacket cannot involve arbitrary 
superpositions of states, being restricted to the space of positive energy solutions. It turns out that a strictly localized 
state cannot be constructed, in contrast with the nonrelativistic case where arbitrarily localized states can be easily 
written down. A relativistic particle is at best "nearly" localized, i.e., concentrated in a small region of space. The 
questions of how well localized the particle can be, and how to characterize these localized states, were dealt with 
in the classical work of Newton and Wigner *9], where a set of natural conditions were stated which any localized 
state should satisfy. The solutions for these conditions are the so-called Newton- Wigner states. These states can be 
characterized as infinite norm eigenstates of a relativistic version of the position operator. In the following, we adopt 
this notion of localizability and use the Newton- Wigner operator to evaluate position averages where necessary. 

The existence of a particular classical limit relies on the existence of an adequate kind of coherent state. So, in 
order to formulate a classical particle limit, one must first look for the analogous of particle- like coherent states in 
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relativistic quantum theory. This problem was addressed in several previous works and we refrain from giving a 
complete list of references here. The states introduced by G. Kaiser in [l(| (which are also particular cases of the 
formalism developed by S. Ali et alii in [2(1]) proved to be a good starting point for our purposes. In order to simplify 
the calculations, we actually work with a simple modification of these states. Before we introduce the coherent states 
we will deal with, let us recall some aspects of the theory of Newton and Wigner. 

Let ip(x, t) be a classical scalar field satisfying the Klein-Gordon equation with mass m in 1 + 1-dimensional 
Minkowski spacetime, and <p(p, uj) its momentum space representation, and assume that <p(p, to) = <j)(p) is restricted 
to the positive mass shell, with uj = \]p 2 + m 2 (we adopt c = 1 throughout the paper). One can view (p(x, t) as a 
first quantized particle wavefunction or as a state in the one-particle sector of the Fock space of the corresponding 
quantum field theory. 

In this context, there are two relevant Hilbert spaces to be considered: Ti 1 = L 2 (R,dp/uj), with the relativistically 
invariant scalar product ((f>\ip) = (4>\^>) n i '■= J R ^ tp{p)ip(p), and Ti 2 — L 2 (R, dp), with the usual scalar product 
(</>|V')'H 2 := f„dp(j>(p)tp(p) (both scalar products interpreted in momentum space representation). The map M^j : 
Ti 2 —>■ Ti 1 defined by (M^j<fi)(p) := t/uj <p(j>) (multiplication operator by \fuj) is unitary. The usual position operator 

on Ti 2 is ifc-§£, and is self-adjoint in some adequate domain. Its counterpart in Ti 1 is q := M ^ (ihj^J the 
so-called Newton- Wigner position operator (at time zero) [T3|. q is also self-adjoint, since M is unitary and 
(q4>)(p) := iti(j^; — 4>{p)- The Newton- Wigner position operator at time t, denoted here by qt, is given by 
q t : = e iut/K qe -iut/K^ for t g Rj j e ; 

d / e -iut/n^\ _ ( d p itp 



and is also self-adjoint in Ti . The momentum operator both in Ti 1 and Ti 1 is just multiplication by p. 

Recall that the Newton- Wigner states @ localized at x at time t, denoted here by ip( x ,t), are given in momentum 
representation in Ti. 1 by {p\tp( x .t)) = , 'l'(x,t){p) = \f^oJ 7 lie l ^ t l h - l P x / h _ They are infinite norm eigenstates of the 
Newton- Wigner operator at time t given in Ij28|) . i.e., qt^(x.t) — xip^ x , t ), what allows to express qt in its spectral 
representation form in Ti 1 as qt — J R x \ipr Xt t)) (V'fei) dx. Notice that J R \i/)( x ,t)) {^P(x.t)\ dx is formally the identity 
operator on Ti. 1 and that (ip( x , t) ^{x 1 , t)) — $( x ~ x ')> as one eas yly checks. The Newton- Wigner wavefunction [l9[ 
associated to a state <j> is given by 

^ w (x, t) := (i> {x , t) \<t>) = _L / ^e-^+^V^m- (29) 



For each t, one has (qt<fi) NW (x, t) — x(f> NW (x, t). Therefore, we can interpret the space of the Newton- Wigner 
wavefunctions <j) NW (x, t) as the spectral representation space of qt, i.e., the space where it acts as a multiplication 
operator [25|. In other words, the wavefunction <p NW (x, t) is the description of the state </> in the spectral representation 
of the Newton- Wigner operator q t . Notice also that {4>\4>') — L 4> NW (x, t) 4>' NW (x, t) dx, and we are allowed to regard 

\(f> NW (x, i)| 2 as the probability density to find a particle at the position x at the instant t. 

With these ingredients we can now define localized versions of the position operator q t associated to (measurable) 
regions O C R by the spectral representation q t ,o '■= Jq x IV^m,*)) {^{x,t) \ dx or, equivalently, by 

(qt,o4>) NW (.x, t) := x X o{x)^ NW {x, t) , (30) 

where Xo{%) is the characteristic function of O. The operator q t: o is self-adjoint and is bounded for bounded O. 

We now turn to the definition of the coherent states introduced by Kaiser in [10j and discuss their more relevant 
properties. For each z := (£ — in, t — ie) € C 2 , with ^, it, t, e 6 IR and e > |7r|, define a coherent state by 

<f> z (p) := N- 1 y/u)exp\i%(T + ie)-i^ + m)] . (31) 
In n J 

The states originally considered by Kaiser differ from those in (f3"Tj) in that the factor ^fuj is absent in his formulation. 
This change corresponds to taking the original states by Kaiser as elements of Ti 2 , with the states in (|3"Tj) being the 
corresponding elements in Ti 1 obtained by applying the unitary map M^. More comments about this are found 
below. The normalization constant N is fixed by the condition 

1 = (4>M = [ ^\MP)\ 2 = 4 / dpe- 2 ^-^' h = ^± Kl (2m\/h), (32) 
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with A = \/e 2 — 7r 2 , where K„, here and below, are the modified Bessel functions of i/-th order (MacDonald's 
functions). This leads to N — 2meKi(2raX / H) / X. Momentum and position averages can be calculated in explicit 
form. 

Their expectation values of p and qt in the states <p z are given by 



<P>*. = 



7T K 2 (2mX/H) 



XK^mX/H) ' 
= € + «(*-t), 



for v — — . 



(33) 
(34) 



From the expressions above we see that the average position of the wavefunction moves with a constant velocity v 
determined by the parameters e, 7r, which also determine the average momentum of the state. The parameter £ is 
the initial position of the coherent state labeled by z, and r the initial instant of time. The usual relativistic relation 
between momentum and velocity is obtained in the K — ► limit, as we will show later. 

It is interesting at this point to discuss the localization properties of the solution of the Klein-Gordon equation 
associated to the coherent states (f> z (p). We will denote these solutions by (p z (x, t). They are given by the Fourier 
transform 



ip g (x, t) = 



'2-K JR2 



dpodpe-^'-^/Vzbo, p)0( Po )S(p 2 - m 2 ) = 



1 



'2-k 



-iujt/h+ipx/h 



4>z(p) 



(35) 



This function depends on x and t only through the combinations x — £ and t—T. Hence, for simplicity, we set r = £ = 
and define zo = (— in, —ie). We are interested in the asymptotic behavior of this wavefunction when x — > oo. In 
order to find it, write 



Vz (x, t) 



dp e 



ipx/h 



1 



-1/2 



/2tt 



exp 



U3 , , V7T 



(36) 



and consider each of the two factors in curly brackets in the integrand independently, so that the transform can be 
computed by the convolution theorem. The Fourier transform of the first factor is just the Newton- Wigner state 
localized at x — 0, 



dp e 



ipx j h 



1 



/2tt 



=w 1/2 oc 



2mh 



1/4 



Ki/4 (-m\x\/h) 



The second factor, after changing to hyperbolic coordinates p = msinhs, to = mcoshs, and substituting 

e + it = pcoshy, n + ix = psmhy, 
with 3?(/o) > and 3(y) G [-7r/2, ir/2], becomes 



rn 
N 



ds exp 



mp 



cosh(s — y) 



cosh s 



2m e + it 
~N~ p 



Ki(—mp/h) . 



Thus, 



f°° e + it ( 1 \ 1/4 

<Pzo(x, t) oc / du K^-mp/h) ( ii ) ^1/4 (-m|x - it|/ft) , 

J— oo P \\ x u \ J 



(37) 



(38) 



(39) 



(40) 



where p is calculated with u replacing x in (|38p . But p ~ u + in for large u, and the MacDonald functions -ftTj, 
display an exponential decay for large arguments. This can be used to prove that the overlap integral in (|40j) has an 
exponential decay for large values of x, i.e., that 



t)\<ne- m W'\ \x\>R, 



(41) 



for some positive n, R. This shows that the coherent state <p z is in fact a wavepacket, localized in some finite region, 
outside of which it falls to zero exponentially with a mass-dependent rate. This asymptotic behavior is characteristic 
of localized relativistic particles. The Newton- Wigner states satisfy exactly the same inequality @. Moreover, the 
notion of essentially localized state in Algebraic Quantum Field Theory is also based on an analog inequality [l[ . 

As we mentioned, the states originally considered by Kaiser differ from those in (|3"Tj) by the factor y/uj. The basic 
properties of the coherent states are not affected by the introduction of this factor, since in both cases the result is 
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a localized wavepacket moving with constant velocity. The advantage in introducing this factor is that the velocity 
v acquires a simple interpretation in terms of the parameters 7r, e, what not only will be useful when dealing with 
the classical limit, but also gives a more direct physical interpretation of these coefficients. Moreover, the expression 
of the associated Newton- Wigner wavefunction is severely simplified, as will be discussed below. The payoff is that 
the spacetime wavefunction ip z (x,t) gets more complicated. The reason for our choice is that in order to define the 
relativistic version of local operators it is natural to use the Newton- Wigner representation, and in consequence of 
this we will work mainly in this representation. 

For the Newton- Wigner wavefunction associated to the coherent states 4> ZQ one gets by the same integral solved 
in (l39l). so that 



1 2m e + it 
N p 



K 1 (mp/h) , 



(42) 



with p given by ([38)1 . The expressions for general £,r are obtained with translations. It follows that the asymptotic 
behavior is given in the new representation by the same expression (|41[) . 

We have proved that the wavefunction decays exponentially with large spatial arguments, for a fixed time. On the 
other hand, from (|34p it is seen that the average position moves with constant velocity. Now we want to discuss how 
the wavepacket spreads about this average motion, so consider the variance of the position distribution. It turns out 
that 



dp 



P -2(t<je-p7r) /1H 



: D{hf 



(43) 



and that limft-^o cr^ = (uniformly for t in compacts). Therefore, the wavepacket is well concentrated about the 
average motion for small h. The momentum is also well-determined in this limit, 



to 4 K 3 {2mX/h) - K x (2m\/h) 
~4 K 1 (2mX/K) 



K 3 (2m\/h) (K 2 {2m\/h) 



Ki{2m\/h) \Ki{2m\/h) 



(44) 



what leads to linifi_>o <Jp = 0. In this limit, p = run/ A = mj(v)v, for 7(f) := (1 — v 2 ) 1 ' 2 , and the usual relativistic 
relation between momentum and velocity is obtained. 



The fact that limn- 
of the coherent states 

iNW 



•o v qt = can be used to determine a nice property of the Newton- Wigner wavefunctions 
Let us now take £ 7^ but r = 0. In the Newton- Wigner representation one has cr? = 

J R dx \<fi z vw (x, i)| z [x — (£ + vt)] 2 . For R > 0, define the time-dependent region Ot, $_ := {x; \x — (£ + vt)\ < R} with 
O c t ^ being its complementary set. Hence, 



R 2 



dx \^ w (x,t) 



< 



dx U N z w {xM[x-{ti + vt)] 2 < ai 



implying that f QC dx \cj)^ w (x,t)\ < D(h) 2 /R 2 . Since lim^^o D(K) — 0, we may claim that, for each fixed R, the 

fraction of the L 2 norm of <f>^ (x, t) outside the time-dependent region O tj ^ is smaller than any desired bound 
Da > as H — > 0. This is a relativistic analog of (TIT)]) . Besides that, it follows by analogous computations that the 
average of the local operator is a good approximation for the average of q t : 



< ^[\Z + vt\+D{H)] 



(45) 



Summing up, the situation is very similar to that found in the nonrelativistic case: coherent states are states 
essentially localized in some region O, outside of which the wavefunction is as small as desired when H — * 0, and 
the average of the position operator can be approximated by a local version. We will now show that the same steps 
followed there can be repeated here, and an approximation scheme can be devised for the free scalar field which leads 
to a classical limit describing a system of relativistic particles. 

First, define iV-body local operators q[ N Q acting on the TV-particle sector of the Fock space generated by Ti. 1 by the 
recipe given in (fl"3)) . At each time t, essentially localized TV-particle states $ w are defined as symmetrized products 
of one-particle coherent states <f>i = (f> Zi situated at disjoint regions Oi — [a, — R, Q>i + R], where a, = £j + i>jt. Let d 
be the smallest distance between these regions, and put A := max{|ai|}. Then, 



1(^1^)1 < 



2D(h) 
R 



(46) 
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and for the operator g tj o k ; 

\mit,oM\ < ( A + R )§^> i{i *j> ( 47 ) 

\(<h,oM < if^fc, (48) 

Kft)^-(ft,o»)* fc | < (A + D(h))^. (49) 
These are the analogs of Eqs. (|T5|l — (p~7|) . Then, it turns out that 

up to an error of order D(h). But the average at the right can be evaluated in the classical limit ft — ► 0, leading to 

This completes our argument. The complete result is the following: there are localized position operators defined 
in each TV-particle sector for the free scalar field whose mean values, when evaluated at essentially localized states 
constructed as the symmetrized products of relativistic coherent states, follow the expected classical trajectories in 
the limit ft — + . 

We have restricted the definition of the local operators to some iV-particle sector, but this is not necessary. One 
can define local operators 0^ =1 qfo in Fock space which measure positions inside a specified region O of space. 
The average of this operator for an essentially localized state \ip) with any number of localization centers Oi is 
approximately zero if Oi H O — for all i, and is approximated by {4>k \ Qt \4>k) if Ok = O for some k. The number 
of particles of the state |i/>} is irrelevant; the local observer at O cannot determine how many particles are there in 
regions of space which are not accessible to it. 

The above results can be viewed in connection with the question of characterizing which quantum field theories can 
be described in terms of particles. The relativistic coherent states we have discussed provide a quantum representation 
for the classical concept of a relativistic particle in free motion. Whenever one uses a quantum field theory to describe 
scattering processes between interacting particles, it is somehow assumed that states which correspond to particles 
do exist. This requirement is usually formulated as the condition of existence of essentially localized states whose 
wavefunction decays exponentially with a mass-dependent coefficient at any fixed instant of time. This is the basic 
idea underlying the approach of [l|. If one is willing to be more restrictive, there is the possibility of imposing the 
additional condition that localized states remain localized at all times, and this line of argument was pursued in 
[I!]. The discussion above suggests an alternative approach. One may try to associate coherent states to particle 
tracks observed in experiments. These states must be essentially localized, but besides that, one should impose the 
condition that the average of suitable local position operators must follow classical trajectories, with a negligible 
dispersion about this average, at least for a time-interval compatible with the experiment in question. In this case the 
question of the existence of a particle interpretation could be recast as the question of existence of a classical limit 
describing particles, and Hepp's analysis would be a natural tool to deal with this problem. 



V. CONCLUSION 



We have defined essentially localized states and local versions of Weyl operators for nonrelativistic bosonic particles, 
and used these concepts to prove the existence of a particle classical limit for a system of N interacting bosonic particles 
in an external potential. The quantum theory of a massive scalar field was proved to have two distinct classical limits, 
describing classical particles or a classical field. Local versions of the Newton- Wigner position operator were defined 
for this theory, and the particle classical limit was obtained as the ft — * limit of the average of these local operators 
evaluated on a class of essentially localized states constructed as symmetrized products of relativistic coherent states 
introduced herein. 

N.Y. thanks FAPESP for financial support. 
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APPENDIX A: TIME-EVOLUTION OF COHERENT STATES FOR LINEARIZED EQUATIONS OF 

MOTION 



A system of N distinguishable particles of the same mass (m = 1, for convenience) is described by the (unsym- 
metrized) Hilbert space Ti,® N of square-integrable functions on R N . Consider the time-evolution generated by the 
time-dependent quadratic Hamilton operator 



N 2 1 

H{t) - E y + 9 E ~ &(*))(* 

t=l 



(Al) 



where the j = 1, . . . , N, describe the classical trajectories of N interacting particles for a solution of the 

corresponding classical system with initial conditions a = [£+i7r]/v2H- In the Heisenberg picture, the time-dependent 
position and momentum operator satisfy linear equations 

ihqi{t) = Pi {t) , m Pi (t) = - E - €s (<))<&(*) . 



identical to the classical equations of motion. These are solved by 

S(t) 



lit) 
Pit) 



where S(t) is a (2N) x (2iV) symplectic matrix (in order that the canonical commutation relations be preserved) 
whose entries depend continuously on t. Writing S(t) in block form as 



-BC 1 



A(t) B{t) 
C(t) D(t) 



1, AB l = BA\CD l = DC 1 ElI. Here, A, B, C, D and 1 are 



the symplectic condition is equivalent to AD* 
N x N matrices, 1 being the identity matrix. 

Now, let the initial state of the system be the coherent state |0) of zero position and momentum in H.® N , and 
ipoix) = (x|0) be the corresponding wavefunction. The time-evolved state at instant t in the Schrodinger representation 



is ^>(x, t) 
equations 



(x| W(t) |0) as usual, where W(t) = Texp -i/h f*dt'H'(t') 



The initial state satisfies the differential 



[ qj +i Pj ]iP (x)=0, j = l,...,N, (A2) 

what in turn implies that [qj(t) + ipj(t)]ip(x,t) = 0, Vj. In terms of the block-components of S(t), one has 

[(A:(t) + iCl(t))x s -h(D r s (t)-iB r s (t))V a ]i>(x,t) = 0. 

The indices r, s = 1, . . . , N label matrix coefficients and the Einstein sum convention is used. Some results about 
symplectic matrices allow for an exact solution of this linear system of differential equations. First, it is advantageous 
to break S into simpler factors, and solve successively for the action of each factor. It was proved in that a 
so-called structured singular value decomposition exists for any symplectic matrix S: one can write 

S = O'DO , (A3) 



where O, O' are orthogonal, D = diag(wi, . . . , ujn, w^f 1 , . . . , w^ 1 ), with ojj > for all j, and the three factors in 
the decomposition are symplectic matrices. That any square matrix M admits a singular value decomposition (i.e., 
M can be written as M — O1D1O2, with 0\ and O2 orthogonal and D\ diagonal with non-negative entries) is a 
well-known result; a structured decomposition is one in each all three factors are symplectic. Second, an orthogonal 
symplectic matrix O has a very simple form, 



= 



U V 

-V u 



(A4) 



where U, V are real matrices such that U — iV is unitary. 
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Hence, for a given time t with |i| < T write S(t) in the form of a singular value decomposition (|A3|) . with O as in 
(|A4[) . Consider first the action of the factor O. One can define new operators 



Y 


= 


<i 


?'. 




P. 



which satisfy the canonical commutation relations [q'j, p' k ] = ihSjk- But, by the Stone-von Neumann Theorem, there 
is, up to unitary equivalence, only one representation of the canonical commutation relations [26]. hence the pair q', p' 
must be unitarily equivalent to the original pair q, p. Let Uo be a unitary operator such that q'j = U^qjUo, and 
P'j = UqPjUo- Use this operator to define the state ip'(x) — Uoip{x) which, from (|A2p . must satisfy 

(Ur-iV s r )(x s + hV s )ip'(x) = 0. 

Multiplying by U 1 + %V t on the left, one gets [x s + HV s )ij, > (x) = 0, whose solution is ip'(x) = ip(x). The initial state 
is not changed by the action of O. The action of the factors D and O' can be studied in the same manner. Define 



Y' 


:= D 


Y 




V" 


:= 0' 


Y~ 


P". 




/_ 




P '" 




p". 



and let Ud, Uo> be unitary operators which accomplish the corresponding unitary transformations. The state ip" = 
Ud^P satisfies the system of differential equations 

(uj r x r + fruj~ 1 X7 r )ij;" (x) = , 

which is solved by the product of Gaussian functions ip"(x) = ^^(Tr/i) -1 / 4 ^^ 2 exp[— (ajjXj) 2 /2h}. This state is not 
changed by the action of O' , so that U o'U v~U o^{x) — ^P"{x). But Uo'UdUo — Us, and the unitary operator which 
correspond to the symplectic matrix S is just the propagator W(t). Therefore, 

i/)(x,t) = J\iTTh)- l/A uj) /2 expl-iujXj) 2 /2h] . 

i 

The time-evolution of a coherent state under a quadratic Hamiltonian can always be solved via a singular value 
decomposition. An algorithm giving the decomposition of an arbitrary symplectic matrix was developed in [22| . 
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